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Abstract 

Let A C Z d be finite. We study the following sandpile model on A. 
The height at any given vertex x € A is a positive real number, and 
additions are uniformly distributed on some interval [a, b] C [0,1]. The 
threshold value is 1; when the height at a given vertex exceeds 1, it 
topples, that is, its height is reduced by 1, and the heights of all its 
neighbours in A increase by l/2d. We first establish that the uniform 
measure [i on the so called 'allowed configurations' is invariant under 
the dynamics. When a < b, we show with coupling ideas that start- 
ing from any initial configuration of heights, the process converges in 
distribution to u, which therefore is the unique invariant measure for 
the process. When a = b, that is, when the addition amount is non- 
random, and a £ Q, it is still the case that [i is the unique invariant 
probability measure, but in this case we use random ergodic theory 
to prove this; this proof proceeds in a very different way. Indeed, the 
coupling approach cannot work in this case since we also show the 
somewhat surprising fact that when a = b £ Q, the process does not 
converge in distribution at all starting from any initial configuration. 



1 Introduction 

The classical BTW- Sandpile Model was introduced by Bak, Tang and Wiesen- 
feld in 1988 as a prototype model for what they called 'self-organized criti- 
cality' [2]. It has attracted attention from both physicists and mathemati- 
cians. In Dhar's paper [4], the model was coined the Abelian Sandpile Model 
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(ASM) because of the Abelian structure of the toppling operators. The 
BTW-sandpile model is defined on a finite subset A of the d- dimensional 
integer lattice. At each site of A, we have a non-negative integer value that 
denotes the height or the number of 'grains ' at that site. At each discrete 
time t, one grain is added to a random site in A; the positions of the ad- 
ditions are independent of each other. When a site has at least 2d grains, 
it is defined unstable and topples, that is, it loses 2d grains, and each of its 
nearest neighbours in A receives one grain. Neighbouring sites that have re- 
ceived a grain, can now become unstable themselves, and also topple. This 
is continued until all heights are at most 2d — 1 again. This will certainly 
happen, since grains at the boundary of the system are lost. See [4, 10, 11] 
for many detailed results about this model. 

As a variant of the BTW model, the Zhang sandpile model was introduced 
by Zhang [13] in 1989. The Zhang sandpile model has some of the flavour 
of the BTW sandpile, albeit with some important and crucial differences: 
(i) in the Zhang-model, the height variables are continuous with values in 
[0, oo) and the threshold is (somewhat arbitrarily) set to 1 (we also speak of 
mass rather than of height sometimes); (ii) at each discrete time, a random 
amount of mass, which is uniformly distributed on an interval [a, b] C [0, 1), 
is added to a randomly chosen site; (iii) when a site has height larger than 1, 
it is defined unstable and topples. In this model this means that it loses all 
its mass and each of its nearest neighbours in A receives a l/2d proportion 
of this mass. 

Since the mass redistributed during a toppling depends on the current 
configuration, the topplings in Zhang's sandpile are not abelian. From simu- 
lations in [13], it is shown that in large volume under stationary measure, the 
height variables of this model concentrate on several discrete 'quasi- values', 
and in this sense, this model behaves similarly to the BTW-sandpile model. 
Recently, this model as well as the infinite volume version of this model have 
been rigorously studied in [5, 6]. 

In the present paper, we discuss a model introduced in [1], which is contin- 
uous like Zhang's sandpile, but abelian like the BTW sandpile. The model 
is related to the deterministic model introduced in [7]. In our model, the 
threshold value of all sites is again set to 1. The only difference with Zhang's 
model, is that when a site topples, it does not lose all its mass, but rather 
only a total of mass 1 instead. Each neighbour in A then receives mass l/2d. 
Note that this is more similar to the BTW sandpile than Zhang's sandpile: 
in the BTW model the number of grains involved in a toppling does not de- 
pend on the actual height of the vertex itself either. In this sense, the newly 
defined model is perhaps a more natural analogue of the BTW model than 
Zhang's model. In this paper we therefore call the new model the CBTW 
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sandpile, where the £ C stands for 'continuous'. 

In this paper, we study existence and uniqueness of the invariant prob- 
ability measure for the CBTW sandpile. In the next section we formally 
define the model, set our notation and state our results. After that we study 
the relation between the CBTW and the BTW-sandpile model, used mainly 
as a tool in the subsequent proofs. 

2 Definitions, notation and main results 

Let A be a finite subset of 7L d and let r\ = (r)(x),x 6 A) be a configuration 
of heights, taking values in X = [0, oo) A . Site x is called stable (in rj) if 
< r)(x) < 1; if 7)(x) > 1, site x is called unstable. A configuration rj is called 
stable if every site in A is stable. We define fl — [0, 1) A to be the collection 
of stable configurations. 

By T x , we denote the toppling operator associated to site x: if A is the 
|A| x |A| matrix 

{1 if x = y, 
otherwise, 

then 

T x rj :=??- A(av). 

Therefore, T x rj is the configuration obtained from r\ by performing one top- 
pling at site x. The operation of T x on rj is said to be legal if rj(x) > 1, 
otherwise the operation is said to be illegal. It is easy to see that toppling 
operators commute, that is, we have 

T x T y rj = T y T x rj. (2.1) 

This abelian property of the topplings implies a number of useful proper- 
ties, the proofs of which are similar to the proof of Theorem 2.1 in [10] and 
are therefore not repeated here. 

Proposition 2.1. Let r] G X and suppose we start the sandpile dynamics 
with initial configuration rj. Then we have that 

1. the system evolves to a stable configuration after finitely many legal 
topplings; 

2. the final configuration is the same for all orders of legal topplings; 
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3. for each iGA, the number of legal topplings at site x is the same for 
all sequences of legal topplings that result in a stable configuration. 



It follows that for rj G X, there is a unique stable configuration rj' G Q 
reachable from rj by a series of legal topplings. We define the stabilization 
operator S as a map from X to Q such that 

Sr] := r]'. (2.2) 

For x G A, u G [0, 1), A u £ denotes the addition operator defined by 

A u xV = S(r ] + u5 x ), (2.3) 

where 

if y ^ x, 
if y = x, 

that is, we add mass u at site x and subsequentially stabilize. According to 
Proposition 2.1, A™ is well-defined, and satisfies 
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A'xA'yi] = AyA^r). 

The evolution of the CBTW model can now concisely be described via 

r H = A%r H - 1 , * = 1,2,3,.-., (2.4) 

where 770 is the initial configuration, X\,X2, ... is a sequence of i.i.d. uni- 
formly distributed random variables on A and Lq, U2, ■ ■ ■ are i.i.d. and uni- 
formly distributed on [a, b], where [a, b] C [0, 1]. The two sequences are also 
independent of each other. 

For 77 G O and W C A, 77] w denotes the restriction of rj to W. For every 
finite subset W of A, r)\w is called a forbidden sub- configuration (FSC) if for 

d\\ x eW 

7}{x) < — :# (nearest neighbours of x in VF), 
2a 

where # denotes cardinality. A configuration rj G ft is said to be allowed if 
it has no FSC's. This notion is parallel to the corresponding notion in the 
BTW-sandpile model — see also below. It is well known that in the BTW 
model, the uniform measure on all allowed configurations is invariant under 
the dynamics (and is the only such probability measure). Let us therefore de- 
note by 1Z the set of all allowed configurations and by fi normalized Lebesgue 
measure on 1Z, that is, the measure which assigns mass such that for every 
measurable B C 1Z we have 
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where Vol(-) denotes Lebesgue measure. 

For every initial probability measure is, z/"' 6 denotes the distribution of the 
process at time t. The first result (similar to the corresponding result in the 
deterministic model in [7]) states not only that fi is invariant for the process, 
but moreover that \x is in fact invariant for each individual transformation A": 

Theorem 2.2. For every x G A,u G [0, 1), A™ is a bijection on 1Z and ji is 
A^-invariant. Hence, /i is invariant for the CBTW-sandpile model. 

When a < b, the situation is as in the traditional BTW model: 

Theorem 2.3. For every 0<a<b<l,/iis the unique invariant probability 
measure of the CBTW model and starting from every measure v on Q, 
converges exponentially fast in total variation to fi as t tends to infinity. 

In the case a = b however, things are more interesting: 

Theorem 2.4. When a = b G [0, 1) but a £ : I = 0, 1, . . . , 2d - 1}, for 

every initial configuration r\ G Q, the distribution of the process at time t 
does not converge weakly at all as t — > oo . 

Theorem 2.5. When a = b and a G : I — 0, 1, . . . , 2d — 1} ; for every 
initial configuration r\ G Q, the distribution of the process at time t converges 
exponentially fast in total variation to a measure pP a as t — > oo. 

In the proof of Theorem 2.5 below, we give an explicit description of the 
limiting distribution \i\ in terms of the uniform measure on the allowed con- 
figurations of the BTW model. In general, the limiting distribution depends 
both on the initial configuration 77 and the value of a. Finally we have 

Theorem 2.6. When a G [0, 1) is irrational and a = b, fi is the unique 
invariant (and ergodic) probability measure for the CBTW model. 

Hence, it is always the case that \x is the unique invariant probability 
measure for the CBTW model, except possibly if a = b G Q. However, 
convergence to this unique stationary measure only takes place when a < b, 
and not (apart from the obvious exceptional case when we start with /i) 
in the case a = b. In the proof of Theorem 2.3, we will use the fact that 
a is strictly smaller than b in order to construct a coupling. The proof of 
Theorem 2.6 is very different from the proof of Theorem 2.3 (although both 
results are about the uniqueness of the invariant measure). In the case a = b 
we use ideas from random ergodic theory rather than a coupling. We could 
have done this also in the proof of Theorem 2.3 but then we would not have 
obtained the exponential convergence corollary. Theorem 2.6 is by far the 
most difficult to prove here, having no coupling approach at our disposal. In 
the next section, we discuss the relation between the CBTW and the BTW 
models. This relation will be used in the subsequent proofs of our results. 
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3 Relation between the CBTW model and 
the BTW model 



We denote by X° and VL° the collection of all height configurations and stable 
configurations in the classical BTW model, respectively. Furthermore, we use 
the notation t x , a x and S° to denote the toppling operator at x, the addition 
operator at x and the stabilization operator in the BTW model, respectively. 

In order to compare the CBTW and BTW models, it turns out to be very 
useful to define the integer and the fractional part of a configuration r\ e X. 

Definition 3.1. For a configuration r] e X , fj is defined by 

rj(x)= [2d V (x)\, (3.1) 
and we call fj the integer part ofrj. We define fj by 

fj(x) = 7](x) mod ^, 

and we call fj the fractional part ofrj. 
The identity 

r, = ±rj + fi (3.2) 

clearly holds. The evolution of the integer part of r\ is closely related to a 
BTW sandpile, while the fractional part is invariant under topplings. These 
two features will make these definitions very useful in the sequel. 

Clearly, r)(x) > 1 if and only if rj(x) > 2d. The operation of taking integer 
parts commutes with legal topplings: 

Lemma 3.2. For r\ e X such that rj(x) > 1, we have 



T X T] = t x T) 

and hence 

Sri = S°fj. 



Proof. For r\ e X, rj(x) > 1 implies r)(x) > 2d. From the toppling rule, we 
have 

T x r] = r]-5 X + 2d S y 

yeA:\y-x\=l 
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which implies that for z G A, 

T^{z) = [2d V (z) - 2d5 x (z) + 5 v^)\- 

yeA:|y— x|=l 

Since both 2dS x (z) and Yl y eA-\y-x\=i $y( z ) are integers, we get 
Tjj(z) = [2dr)(z)\ - 2d5 x (z) + £ 5 y (z) 

yeA:\y-x\=l 

which is equal to (t x rj)(z). Therefore we obtain T x i] = t x fj. It now also follows 
that Srj = S°rj. □ 

For a configuration £ G Af , let C(£) be the set 

C(£) := {77 G A- : 77 = 6}, 

that is, C(£) is the set of all configurations which have £ as their integer part. 
The sets C(£), £ G partition f2. Indeed it is easy to see that for £, £' G f2° 
with £ 7^ £', we have 

C(£) n C(£') = (3.3) 

and 

n=yc(0- (3-4) 

Before continuing to the next observation, we recall the definition of a 
BTW- forbidden sub-configuration. For £ G f2°, £|vk is a FSC if for all x G W 

£(x) < #( nearest neighbours of x in W). 

By 72.°, we denote the set of allowed configurations (that is, configurations 
without FSC's) in the BTW model, which is also the set of all recurrent 
configurations, see [10]. For rj and its corresponding integer part 77, we have 
the following simple observation, the proof of which we leave to the reader. 

Lemma 3.3. For rj G Q, rj is allowed if and only ifrj is allowed in BTW. 

Definition 3.4. A configuration r] G Q is called reachable if there exists a 
configuration rf G X with rj'(x) > 1 for all x G A, such that 

rj = Srj' . 

We denote the set of all reachable configurations by VJ . 
Theorem 3.5. K = K'. 
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Proof. Let rj G 1Z. By Lemma 3.3, rj is BTW-allowed and therefore also 
recurrent (Theorem 5.4 in [10]). From [10, Theorem 4.1], we have that for 
every x G A, there exists n x > 1 such that 

a^rj = rj. 

Therefore 

H a 2 x dn ^ = S°(rj + J2 2dn ^) = % (3-5) 

Now let 

be a configuration in X°; note that £(x) > 2d for all x G A. Let 

»/ = ^ + * (3-6) 

Then 77' G X and 

1. r/(:r) > 1 for all x G A; 

2- ? = 

3. 77' = rj. 

We will now argue that Srj' = rj. /From Lemma 3.2 and (3.5) we have 

Srf = S°rf = S°£ = rj. 
Furthermore, since fractional parts are invariant under S, we have 

Srf = r]' = fj, 

and hence Srj' = rj. 

For the other direction, if 77 G 1Z', there exists an rj' G X with rj'(x) > 1 
for all x G A, and such that 

rj = Srj'. (3.7) 
Clearly, r)'(x) > 2d. By Lemma 3.2, we can now write 

xeA 

where ^ max (x) = 2d — I, for all igA. This means that Srj' is a configuration 
obtained by additions to £ max . Since £ max is allowed for the BTW model, by 
Theorem 5.4 of [10], ^ max is also recurrent. But then Srj' is also recurrent 
and therefore allowed. It then follows that Srj' is allowed for the CBTW 
model. □ 
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Corollary 3.6. For x G A, w G [0, 1), 1Z is closed under the operation of ' A u x . 

Proof. For rj G TZ, by Theorem 3.5, there is a rj' G A" with 7/(2;) > 1 for all 
iGA, such that 

77 = 1S7/. 

By the abelian property of toppling operators, 

A u x rj = S(Srj' + u5 x ) = S(rj' + u5 x ) 

and 77' + is a configuration with only unstable sites. It follows that 
A u x r] G K', and hence A u x r) eK. □ 

Lemma 3.7. (1) Vol(ft) = det(A); 
^ For every f G //(C(0) = ^ 

Proof. (1) We have 

Vol(tt) = ^ Vol(C(0). 

For each £ G ft , Vol(C(f)) = (2d)-l A l, hence Vol(ft) = |ft°|(2d)-l A l In the 
BTW-sandpile model, it is well known that \TZ°\ = det(A°), where A° is the 
toppling matrix, see [10, Theorem 4.3]. Since A° = 2dA, we have 

Vol(ft) = (2rf)l A l det(A)(2rf)~l A l = det(A). 

(2) is immediate from the definitions. □ 

4 Proof of Theorem 2.2 

We denote the sites in A by xi, x 2 , ■ ■ . , X|a| and define the collection £ as 
£={[c,d)= J] [c k ,d k ):0<c k <d k <^Y 

l<fc<|A| 

where c = (ci, . . . , cja|) and d = (d±, . . . , d|A|)- Note that £ is a 7r-system. 
For £ G ft , let 

C(£, [c,d)) = { v eC(0--ve[c,d)}, 

which is the set of configurations whose integer part is £ and whose fractional 
part is in the interval [c, d). Let 

I={C^,[c,d)):[c,d)eC^en°}. (4.1) 
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X is also a 7r-system. In order to show that fi is ^"-invariant, it suffices to 
show that 

Mr? : A u xV EB} = ^B) (4.2) 

for all B G X, and we will do that by direct calculation. 

Let B = C(£, [c, d)) and ( G -B, and define the configuration r\ by 

fj(z) = (C(z) - u5 x (z)) mod ^, (4.3) 

and 

f (a -i )LM »j e if C »> Mmo d^; 

\(aJ 1 ) L2dttJ+1 ^ ^ C(z) < « mod l ' J 

We claim that this rj is the unique rj e 7Z such that A"?? = £. To see this, 
first note that for any r\ G 71, fj differs from A^r] only at the site x, and that 

A%r)(x) = (t](x) + u) mod — . 

2d 

This shows that for any r] G 7Z such that A^rj = (, its fractional part is given 
by (4.3). Now there are two possibilities: 

1. (A^r])(x) = ((x) > u mod 

2. (A^)(rr) =C(x) <wmod^. 

A little algebra reveals that in the first case, 

A^ = S°(r]+[2du\5 x )=ai 2du ^, 

and in the second case, 

~A%q = S°{fj+([2du\ + 1)4) = a^ +1 rj. 

In words, depending on the fractional part of ( at site x, the addition of u 
to site x corresponds to either adding \2du\ or \2du\ + 1 particles at x in 
the BTW model. It follows that any r) G 71 such that A^r/ = ( must have 
integer part given by (4.4). We conclude that for every (6 5, the 77 defined 
by (4.3) and (4.4) is the unique 77 G TZ such that A"?7 = (. It follows that A™ 
is a bijection. 

Next we show that A™ preserves the measure /i. For this, note that 
the inverse image of B naturally partitions into two sets: if x — Xi, define 
[c 1 , d 1 ), [c 2 , d 2 ) to be intervals that differ from [c, d) in the i th coordinate only, 
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to the effect that [c},dj) = [c i: di) fl [u mod ^,^) and [cf,df) = [cj,dj) fl 
[0, u mod Note that [c, <f) = [c 1 , (i 1 ) U [c 2 , rf 2 ). From the above we have 

: A> G C(£, [c\d 1 ))} = C((a~T dui Z, [c',d')) 

and 

{ V : A> G C(£, \c\d 2 ))} = C((a?)W +1 t, [c",d")), 

where [cf, d') and [c", d") are intervals that differ from [c, d) only in the 
jth coordinate to the effect that [c'^d 1 ^ = [c\ — u mod j^,dj — u mod ^) 
and [<%, d'l) = [c 2 i +±-u mod ±^ + ±- u mod i). Finally, 

Vol{C((a; 1 ) [c', } = Vol{C(£, [c 1 , rf 1 )) } 

and 

Vo\{C((a-T du ^, [c'\d"))} = Vol{C(£, [c 2 ,rf 2 ))}, 
which implies that ji is ^"-invariant. □ 



5 Proof of Theorem 2.3 

We prove Theorem 2.3 with a coupling. Let ?), ( G be two initial con- 
figurations and let rj t ,(t be two copies of the CBTW sandpile starting from 
r], ( respectively (hence r\ = r] and ( = Co)- Addition amounts at time t 
are Uf, respectively , the addition sites are X^,X^ respectively. All these 
random quantities are independent of each other. In the proof, we will couple 
the two processes, and in the coupling, random variables will be written in 
'hat'-notation; see below. 

For every x G A and t — 0, 1, . . ., let D t (x) be defined as 

A(x) = T 4- T (ifc(x)-C«(x)). (5.1) 

I b-a I 

We now define a coupling of the CBTW realisations starting from 77 and ( 
respectively. The sites to which we add are copied from the 77-process, that 
is, we define 

X? = X< = X? (5.2) 

and write X t (without any superscript) for the common value. The addition 
amount in the ^-process is unchanged, that is, 

V? = V?, 
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but in the (^-process we define 

U} = [U? + D (X t ) - a] mod (b - a) + a. 

Hence, in the coupling the ^-process evolves as in the original version, but 
the ^-process does not. It is not hard to see that this definition gives the 
correct marginals and that if U? G [^^, ^p], then 

U} = U? + -J-fopT*) - Co(X t )). (5.3) 

I b-a I 

Let us now say that event O occurs if 

1. between times 1 and 1 + |A||~^] (inclusive), all sites are chosen as 
addition sites exactly times; 

2. between times 1 and 1 + |A| (inclusive), the addition amounts in 
the 77-process are all contained in the interval [^^, 9 ^-]- 

If O occurs, then we claim that at time 1 + |A| [^], the two processes 
are in the same state. Indeed, by the abelian property we can obtain the 
configuration at time t by first adding all additions up to time t, and after that 
topple all unstable sites in any order. If O occurs, and we defer toppling to 
the very end, all heights will be the same in the two processes by construction, 
and hence after toppling they remain the same. 

The probability that O occurs is uniformly bounded below, that is, uni- 
formly in the initial configurations rj and (. Indeed, all that is necessary is 
that all sites are equally often addition sites, and that the addition amounts 
are in the correct subinterval of [a, b]. These two requirements are indepen- 
dent of the starting configurations. Of course, the evolution of the coupling 
does depend on the initial configurations via the relations between the added 
amounts. That is where the subtlety of the present coupling lies. 

Hence, to finish our construction, we first see whether or not O occurs. If 
it does we are done. If it does not, we start all over again, with the current 
configurations at time t as our new initial configurations, and D t (x) instead of 
D (x). Continuing this way, we have a fixed positive probability for success 
at each trial and therefore the two processes will almost surely be equal 
eventually. Due to the fact that the success probability is uniformly bounded 
below, convergence will be exponentially fast. The result now follows by 
choosing i] and ( according to the distributions v and /i, respectively. □ 
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6 Proof of Theorems 2.4 and 2.5 



We start with the proof of Theorem 2.4. 
Proof of Theorem 2.4- Define the function g by 



g(rj) = expUdni^2r]( 



,{x 



It is easy to see that g is continuous and, of course, bounded. Denoting the 
distribution of the process at time t by v t , for v t to converge weakly to v, 
say, it must be the case that 

J gdv t -> J gdv. (6.1) 

However, since at each iteration of the process, exactly one of the fractional 
parts is increased by a and subsequently taken modulo l/2d, we have that 
Yl,x£k r lt{ x ) m °d l/2d is //(-surely equal to (XLeA Vo( x ) + ta) mod l/2d, and 
hence the sequence of integrals on the left of (6.1) does not converge at all, 
unless a is a multiple of l/2d. □ 

We now move on to the proof of Theorem 2.5. In order to describe 
the limiting measure //JJ appearing in the statement of the theorem, we first 
introduce the following notation. For a measure v on the space of height 
configurations X, Sv denotes the measure on the stable configurations Q 
defined by 

Sis(B) = v{7] E X : Sr] E B}, 
for every measurable set B C Q. 

Proof of Theorem 2.5. When a — 0, rjt — r] for all t = 1,2, ... , and the 
limiting distribution /Iq is point mass at rj. 

So assume I > and take a = l/2d. We first consider the case in which 
r){x) =0 for all sites x, and proceed by a coupling between the BTW and 
the CBTW model. To introduce this coupling, let (X t ) denote a common 
sequence of (random) addition sites, and define 

t t 

s=l s=l 

This couples the vector Q t of all additions until time t in the CBTW model, 
with the vector Q° of all additions until time t in the BTW model. 
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In order to arrive at the BTW configuration at time t, now we can first 
apply all additions, and then topple unstable sites as long as there are any. 
In doing so, we couple the topplings in the BTW model with those in the 
CBTW model: if we topple a site x in the BTW model, we topple the 
corresponding site in the CBTW model / = Ida times. If we topple at x, 
then in the BTW model, x loses 2d particles, while all neighbours receive 1. 
In the CBTW model, x loses total mass Ida while all its neighbours receive a. 
Hence, the dynamics in the CBTW model is exactly the same as in the BTW 
model, but multiplied by a factor of a. In particular, since all the topplings 
in the BTW model were legal, all the topplings in the CBTW model must 
have been legal as well. Furthermore, when the BTW model has reached the 
stable configuration n t = S°0°, then the corresponding configuration in the 
CBTW model is simply an t . However, an t need not be stable in the CBTW 
model. Hence, in order to reach the CBTW configuration at time t, we have 
to stabilize an t , leading to the stable configuration p t = S9 t . 

From the results in [10], we know that the distribution of K t converges 
exponentially fast (as t — > oo) in total variation to the uniform distribution 
on 1Z°. Hence, pt is exponentially close to Su a , where v a is the uniform 
measure on the set {a£ : £ G 1Z°}. 

This settles the limiting measure in the case where we start with the 
empty configuration. If we start with configuration r\ in the CBTW model, 
we may (by abeliannes) first start with the empty configuration as above- 
coupled to the BTW model in the same way — and then add the 'extra' rj 
to p t at the end, and stabilize the obtained configuration. It follows that the 
height distribution in the CBTW model converges to Sis^, where v% is the 
uniform measure on the set {rj + a£ : £ G 1Z°}. □ 

7 Proof of Theorem 2.6 

The proof of Theorem 2.6 is the most involved. It turns out that the view- 
point of random ergodic theory is very useful here, and we start by reformu- 
lating the sandpile in this framework; see [8] for a review of this subject. 

We consider the CBTW with a = b, that is, with non-random additions at 
a randomly chosen site. As before, we denote by p the uniform measure on 71. 
For every x G A, we have a transformation A a x : 1Z — > 1Z which we denote 
in this section by A x (since a is fixed). Recall that each A x is a bijection 
by Theorem 2.2. The system evolves by each time picking one of the A x 
uniformly at random (among all A x , x G A) independently of each other. 
The product measure governing the choice of the subsequent transformations 
is denoted by p, that is, p assigns probability 1/|A| to each transformation, 
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and makes sure that transformations are chosen independently. The system 
has randomness in two ways: an initial distribution v on TZ and the choice 
of the transformations. To account for this we sometimes work with the 
product measure v x p. 

A probability measure v on TZ is called invariant if 

that is, v preserves measure on average, not necessarily for each transforma- 
tion individually. We call a bounded function g on TZ v-invariant if 

' ' rrGA 

for //-almost all rj E TZ. We call a (measurable) subset B of TZ invariant 
if its indicator function is a zz-almost everywhere invariant function; this 
boils down to the requirement that up to sets of //-measure 0, B is invari- 
ant under each of the transformations A x individually. Finally, we call an 
invariant probability measure v on TZ ergodic if any //-invariant function is 
a z/-a.s. constant. These definitions extend the usual definitions in ordinary 
(non-random) ergodic theory. It is well known and not hard to show (see 
[8, Lemma 2.4]) that v is ergodic if and only if every invariant set has im- 
probability zero or one. As before, we denote by fj and f] the integer and 
fractional parts of the configuration r\. 

Lemma 7.1. Let A be an invariant probability measure, B a \-invariant set 
and C a set such that X(C A B) = 0. Then ^(A^C A B) = for all x e A. 

Proof. From the assumptions and invariance of A, it follows that the sets 
A~ X B A B, A~ 1 (B \ C) and A^C \ B) all have A-measure 0. Now suppose 
cu E A~ l C \ B. Then either u E A~ X B \ B or else u E A-\C \ B). Since 
both these sets have A-measure 0, A(A a T 1 C \ B) — 0. Next suppose that 
u) E B \ A~ l C. Then either u E B \ A~ X B or else u E A~ 1 {B \ C). Again, 
since both these sets have A-measure 0, \(B \ A~ X C) = 0. □ 

A version of the following lemma is well known in ordinary ergodic theory, 
see e.g. [3, Proposition 5.4]. 

Lemma 7.2. Let v be an ergodic probability measure (in our sense) and 
let A be an invariant probability measure which is absolutely continuous with 
respect to v . Then A is ergodic, and therefore A = v . 
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Proof. For u G Z> , write 

T u :=\{Af x \ 

xeA 

By commutativity, the order of the composition is irrelevant. Now suppose 
that B is a A-invariant set. Define 



C == n U Tu l B- 



n=O ueZ A n 

In words, C is the set of configurations r\ such that for all n, there is a w with 
u(x) > n at every ieA, for which T u i] G -B. Since T u (A x n) = T u+ s x r] ) it is 
not difficult to see that A x r\ G C if and only if i] G C. Hence A~ X C = C for 
all x G A. In particular, C is a //-invariant set, so by ergodicity of v>, either 
v{C) = or v(C c ) = 0. But since A is absolutely continuous with respect 
to u, this implies that either A(C) = or A(C C ) = 0. 

To prove that A is ergodic, it therefore suffices to show that X(B) = A(C), 
or equivalently, \(C A B) =0. To this end, we first claim that 

(CAB)c |J (T~ X B A B). (7.1) 

Indeed, if r\ G C \ 5, there must be a w G Z> such that T u i] G -B, hence 
r] G T^ 1 !? \ B. And if rj G B\C, then there are only finitely many n such 
that T M r] G 5 with = n for all x G A. Hence rj E B \ C implies that 
for some n > we have T u r] (jL B with u[x) = n for all x G A, and therefore 
i] E B \ T~ 1 B for this particular w. This establishes (7.1). 

To complete the proof, note that by repeated application of Lemma 7.1, 
it follows that for every u G Z£„, X^^B A B) = 0. Hence A(C A B) = 
by (7.1), and we conclude that A is ergodic. Now it follows from Theorem 2.1 
in [8] that v x p and A x p are also ergodic (in the deterministic sense w.r.t. 
the skew product transformation). But since A is absolutely continuous with 
respect to u, A x p is absolutely continuous with respect to v x p, and hence 
must be equal to v x p by ordinary (non-random) ergodic theory. It follows 
that v — A. □ 

We will apply Lemma 7.2 with the uniform measure fx in the role of v. Be- 
fore we can do so, we must first show that fx is ergodic, which is an interesting 
result in its own right. 

Theorem 7.3. When a = b Q, the uniform measure fx on 1Z is ergodic. 
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Proof. We denote the sites in A by x±, x 2 , ■ ■ ■ , x m , and identify a configuration 
i] G [0, oo) A with the point (r)(xi), . . . , r)(x m )) in [0, oo) m . By Proposition 3.1 
in [10], there exists a n = (rii, . . . , n m ) G Z"\ such that for alH = 1, 2, . . . , m, 



Now let A be the rectangle 



Wo' 



• X 



0, 



2d 



and denote Lebesgue measure on A by A. Write Tj for the translation by a 
modulo rii/2d in the i* h coordinate direction on *4. 

A subset D of „4 is called \-invariant if Xij^D) = \(D) for all Tj. We 
claim that for any A-invariant set D, either X(D) = or X(D C ) = 0. Although 
this fact is probably well known, we give a proof for completeness. Write Id 
for the indicator function of D, and for k G Z m , denote by c k the Fourier 
coefficients of Id- Then 

c k = / I D (u)fk(u)d\(uj) = / f k (u)d\(uj), 

J A JD 



where 



2m — k i LU i 



u e A. 



Now suppose that fcj ^ for some 1 < j ' < m. Then, changing coordinates 
by applying Tj, we have that 

f f 2ni — k a 4d7ri—a 

Ck= fk{rjuj) d\(Tju) = / fk(v)e 3 dX(u) = c k e n i , 

JtT^D Jd 

because Lebesgue measure is invariant under Tj and D is A-invariant. But 
since a is irrational, this implies that = 0. It follows that Id = c Q A-almost 
everywhere. Hence c G is either or 1, so that either X(D) = or X(D C ) = 0. 
Next note that A is composed of the cubes 



Ck 



2d' 2d 



x 



k 2 k 2 



2d' 2d 



x • • • x 



k k 



m n m 



2d 2d 



where k G Z m with < ki < rij for % — 1,2, ... ,m. Using this fact, we define 
a map tp : A — > 1Z as follows. If a; is in the cube C^, let 



&=n™ 1 °s^ ,M ' 
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where £ max (x) = 2d - 1 for all x e A, and set 

Thus, if) simply translates the cube Ck onto the cube of configurations in TZ 
whose integer part is Notice that multiple cubes in A may be mapped 
by ip onto the same cube in TZ, but that -0 : ^4. — > 72. is surjective, because 
every allowed configuration of the BTW model can be reached from £ max 
after a finite number of additions (and subsequent topplings). Furthermore, 
it is easy to see that a translation on A corresponds to an addition A x . 
on 1Z, in the sense that 

^{nuj) = A x ^(lu). (7.2) 

Now suppose that B C 1Z is /i-invariant. Since \i is normalized Lebesgue 
measure on 7Z and A is Lebesgue measure on A, we have that 

xty-^B)) = K^\B) n c k ) = Voi(K)J2v(B n ^(c k )), 

k k 
where the sum is over all cubes in A. Because n(A~^B A B) = 0, this gives 

\{r\B)) = voi(rc) 5>(4^ n ^(c fe )) = £ A(r Ha- 1 ^) n c fc ). 

By (7.2), we have that 

and it follows that A^ -1 ^)) = A(r. _1 ^~ 1 (i?)), hence is a A-invariant 
set. Therefore, either A(f/> _1 (-B)) = or A(f/> _1 (.B C )) = 0. From the construc- 
tion it then follows that either /i(B) = or fi(B c ) = 0. Therefore, ji is 
ergodic. □ 

The next lemma will be used to deal with the evolution of the joint 
distribution of the fractional parts. In order to state it we need a few defi- 
nitions. Consider the unit cube I m = [0, l] m . Let a be an irrational number 
and let 9 be the measure which assigns mass 1/m to each of the points 
(a, 0, . . . , 0), (0, a, 0, . . . , 0), . . . , (0, . . . , 0, a). Denote by 9* n the n-fold convo- 
lution of 9, where additions are modulo 1. Translation over the vector x 
(modulo 1 also) is denoted by r x . We define the measure \x x N by 



1 N ~ x 
x — __ 9* n r~ 1 
N ^ 



n=0 
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In words, this measure corresponds to choosing n G {0, 1, . . . , iV — 1} uni- 
formly, and then applying the convolution of n independently chosen trans- 
formations with starting point x. 

Lemma 7.4. For all x, ijl x n converges weakly to Lebesgue measure on I m as 
N -> oo. 

Proof. The setting is ideal for Fourier analysis. It suffices to prove (using 
Stone- Weierstrass or otherwise) that 

J fkdfx x N -> J f k d\, 

where A denotes Lebesgue measure and 

f k (y) = e 2 ™ k -y = e w ST=i'i«i, 

for k = (k±, . . . , km) G Z m and y G I m . 

When k = (0, . . . , 0), all integrals are equal to 1. For k ^ (0, . . . , 0), 
/ f k d\ — 0, and hence it suffices to prove that for these k, 

J f k dfi x N -> 0. 

For n > and I = (/i, . . . ,l m ) G Z> such that /i + \- l m = n, 

9* n r- 1 {x + la} = —( ' , n A 
Hence, by the multinomial theorem, 

^ i 1+ ...+i m =n m Vl.«2,...,W J=1 

(m , \ n 

E- e27r ^J =:/*(*)(«*)"■ 

Note that < 1 and a k 7^ 1 because a is irrational. Therefore, 

/1 y - . 1 11 z' w 

A « = /.Wtt J>*) b = AOOaF i_* "> °' 
iV n=0 " fc 

as required. □ 
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We will use Lemma 7.4 to understand the fractional parts in the sand- 
pile, working on the cube [0, l/2d] m for m = |A| rather than on [0, l] m . Since 
topplings have no effect on the fractional parts of the heights, for the frac- 
tional parts it suffices to study the additions, and a point x G [0, l/2d] m 
corresponds to all fractional parts of the m sites in the system. 

Proof of Theorem 2.6. We have from Theorem 7.3 that the uniform mea- 
sure /i on 1Z is ergodic. Suppose now that u is another invariant measure. 
We will show that ^{B) = implies v(B) = 0, and according to Lemma 7.2 
it then follows that v = \i. 

Consider the "factor map" / defined via f(rj) = fj, that is, / produces 
the fractional part when applied to a configuration. The map / commutes 
with the random transformations we apply to r], in the sense that 

f(Mv)) = Mfiv)), 

where A x (fj) is the configuration of fractional parts that results upon adding a 
to the height at site x. Write v = vf^ 1 and jl = /i/ -1 . Because / commutes 
with the random transformations and v is invariant, for any measurable 
subset B of [0, l/2rf]l A l we have that 

<r\B)) = ^Y, < A * l r\B)) = JL £ v(r i (X- i B)). 



Therefore, we have 

' ' xeA 

and hence v is invariant. We claim that \x — v. To prove this, it suffices to 
show that fl(R) = v(R) for any rectangle R. Since rectangles are /i-continuity 
sets, by Lemma 7.4 and bounded convergence, 

p,(R) — J fi(R)dv = J hm J I R d(j,% dv(h) = hm J J I R dfj,%du(h). 

Here 11% is the analogue of the measure defined in Lemma 7.4 on the space 
[0, l/2d]' A ' rather than [0, l] m ; it corresponds to choosing n G {0, 1, ... , A^— 1} 
uniformly, and then applying n uniformly and independently chosen trans- 
formations to the fractional height configuration h. Thus we also have that 

J j i R dfi h N du(h) = J 1 ^ 7 «(A« • ■■ A v» h ) d *( h ) 

1 iV " 1 1 

n=0 y 1 ,...,y n eA 
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By invariance of v, this last expression is equal to 

N-l 
n=0 

and we conclude that jl(R) = v{R) for every rectangle R, hence jl — v. 

For the final step, let B be such that n(B) = and write B = f(B). We 
claim that then jl(B) = 0. Indeed, the inverse image of B is a collection of 
points of the form x + (2d)~ 1 r], where x G B and r] G Z A . If /i(B) = 0, then 
also this collection of points has /x-measure 0. Since jl — i>, it follows that 
v{B) = 0, hence also v{B) = 0. This concludes the proof. □ 

Acknowledgment. It is a pleasure to thank Michael Keane for very inter- 
esting discussions. 
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